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Abstract. We report on a mechanism which may lead to 
a spin-up of the surface of a rotating single star leaving 
the Hayashi line, which is much stronger than the spin-up 
expected from the mere contraction of the star. 

By analyzing rigidly rotating, convective stellar en- 
velopes, we qualitatively work out the mechanism through 
which these envelopes may be spun up or down by mass 
loss through their lower or upper boundary, respectively. 
We hnd that the first case describes the situation in re- 
treating convective envelopes, which tend to retain most 
of the angular momentum while becoming less massive, 
thereby increasing the specific angular momentum in the 
convection zone and thus in the layers close to the stellar 
surface. We explore the spin-up mechanism quantitatively 
in a stellar evolution calculation of a rotating 12 M Q star, 
which is found to be spun up to critical rotation after 
leaving the red supergiant branch. 

We discuss implications of this spin-up for the circum- 
stellar matter around several types of stars, i.e., post-AGB 
stars, B[e] stars, pre-main sequence stars, and, in partic- 
ular, the progenitor of Supernova 1987A. 

Key words: stars: evolution - stars: rotation - circum- 
stellar matter - supernova 1987A 



1. Introduction 

The circumstellar matter around many stars shows a re- 
markable axial symmetry. Famous examples comprise Su- 
pernova 1987A (Plait et al. 1995; Burrows et al. 1995), 
the Homunculus nebula around 77 Carina and other nebu- 
lae around so called Luminous Blue Variables (Nota et al. 

1995) , and many planetary nebulae (Schwarz et al. 1992). 
A less spectacular example are B[e] stars, blue supergiants 
showing properties which might be well explained by a cir- 
cumstellar disk (Gummersbach et al. 1995; Zickgraf et al. 

1996) . Many of these axisymmetric structures have been 
exnlained in terms of interacting' winds of rotatinsr stars 



(cf. Martin & Arnett 1995; Langer et al. 1998; Garcia- 
Segura et al. 1998), which may be axisymmetric when the 
stars rotate with a considerable fraction of the break-up 
rate (Ignace et al. 1996; Owocki et al. 1996). However, up 
to now only little information is available about the evo- 
lution of the surface rotational velocity of stars with time, 
in particular for their post-main sequence phases. 

Single stars which evolve into red giants or supergiants 
may be subject to a significant spin-down (Endal & Sofia 
1979; Pinsonneault et al. 1991). Their radius increases 
strongly, and if the specific angular momentum were con- 
served in their surface layers (which may not be the case; 
see below) they would not only spin down but they would 
also evolve further away from critical rotation. Moreover, 
they may lose angular momentum through a stellar wind. 
Therefore, it may appear doubtful at first whether post- 
red giant or supergiant single stars can retain enough an- 
gular momentum to produce aspherical winds due to ro- 
tation. 

However, by investigating the evolution of rotating 
massive single stars, we found that red supergiants, when 
they evolve off the Hayashi line toward the blue part of 
the Hertzsprung-Russell (HR) diagram may spin up dra- 
matically, much stronger than expected from local angular 
momentum conservation. In the next Section, we describe 
the spin-up mechanism and its critical ingredients. In Sec- 
tion 3 we present the results of evolutionary calculations 
for a rotating 12 M Q star, which provides a quantitative 
example for the spin-up. In Section 4 we discuss the rel- 
evance of our results for various types of stars, and we 
present our conclusions in Section 5. 

2. The spin-up mechanism 

2.1. Assumptions 

In this paper, we discuss the evolution of the rotation rate 
of post-main sequence stars, i.e. of stars which possess 
an active hydrogen burning shell source. This shell source 
converts hydrogen into helium and thus increases the mass 
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the star, and for the following argument we can consider 
the shell source as fixed in Lagrangian mass coordinate. 

The shell source provides an entropy barrier, which 
separates the high entropy envelope from the low entropy 
core, and it marks the location of a strong chemical discon- 
tinuity, i.e. a place of a strong mean molecular weight gra- 
dient. Both, the entropy and the mean molecular weight 
gradient, act to strongly suppress any kind of mixing 
through the hydrogen burning shell source. This concerns 
chemical mixing as well as the transport of angular mo- 
mentum (e.g., Zahn 1974; Langer et al. 1983; Meynet & 
Maeder 1997). Therefore, in the following, we shall regard 
the angular momentum evolution of the hydrogen-rich en- 
velopes of the stars under consideration as independent of 
the core evolution. 

This picture is somewhat simplified, as due to the inhi- 
bition of angular momentum transport through the shell 
source the post-main sequence core contraction and enve- 
lope expansion results in a large gradient in the angular 
velocity at the location of the shell source, i.e. a large shear 
which may limit the inhibiting effects of the entropy and 
mean molecular weight gradients. However, with the cur- 
rent formulation of shear mixing in our stellar evolution 
code, we find the angular momentum transport to be in- 
significant (cf. Sect. ||). In any case, the total amount of 
angular momentum in the helium core is much smaller 
than that in the hydrogen-rich envelope ( <J 2 % in the 
12 M© model discussed below), so that even if all of that 
would be transported out into the envelope it would not 
alter its angular momentum balance much. 

As a star evolves into a red supergiant, its envelope 
structure changes from radiative to convective, starting at 
the surface. If we assume the envelope to be in solid body 
rotation initially, it can not remain in this state without 
any transport of angular momentum, unless 
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remains constant with time everywhere, i.e. r 3 p remains 
constant. Here, m is the Lagrangian mass coordinate, r 
the radius, p the density, and i = k 2 r 2 is the specific 
moment of inertia. For a sphere of radius r the gyration 
constant is k 2 = 2/3, but even for deformed equipotential 
surfaces of average radius r one still finds k 2 to be of order 
unity. The assumption that k 2 does not depend on r was 
used to derive Eq. ([!]). Obviously, the above homology 
condition does not hold for stars which change their mode 
of energy transport in the envelope, since the polytropic 
index n changes between n ~ 3 in the radiative envelope 
and n ~ 1.5 in the convective state. 

In the following we assume that convection tends to 
smooth out angular velocity gradients rather than angu- 
lar momentum gradients, i.e. that convective regions tend 
to be rigidly rotating. This is certainly a good approxima- 



time scale (r conv := Hp /v conv ), and it may also hold for 
more rapid rotation if convective blobs can be assumed 
to scatter elastically (cf. Kumar et al. 1995). Latitudi- 
nal variations of the rotational velocity as deduced from 
hclioscismological data for the solar convective envelope 
(Thompson et al. 1996) are not taken into account in our 
ID stellar evolution calculation; however, the latitudinal 
averaged rotation rate of the solar convection zone devi- 
ates by no more than 5 % from solid body rotation (cf. 
e.g. Antia et al. 1997). 



Although, for the considerations in Sect. 2.2 we assume 



rigid rotation to persist in convection zones, the neces- 
sary condition to make the spin-up mechanism described 
here work is only that convection transports angular mo- 
mentum on a time scale which is short compared to the 
evolutionary time scale of the star, and that it leads to a 
characteristic angular momentum distribution in between 
the cases of constant angular velocity and constant angu- 
lar momentum. The efficiency of the spin-up is largest for 
constant angular velocity and drops to zero for the case of 
constant angular momentum. The mechanism we present 
here is not restricted to convection, but any transport of 
angular momentum that acts towards solid body rotation 
is suitable to accomplish what we describe here. 

2.2. The spin-down of convective envelopes 

The rotation frequency of a rigidly rotating convective en- 
velope depends on its moment of inertia and its angular 
momentum. Both are altered by mass loss from this en- 
velope, the former by loss of mass, the latter by the ac- 
companied loss of angular momentum. Additionally, the 
envelope's moment of inertia is also affected by changes of 
its density stratification. 

Here, we want to discuss two processes which can 
change the rotation frequency of a rigidly rotating con- 
vective envelope without employing global contraction or 
expansion of the star, but rather by mass outflow through 
its upper or lower boundary. We will show that the first 
case leads to a spin-down, while the latter spins the enve- 
lope up. 

The spin-down of mass losing rigidly rotating en- 
velopes can be understood by breaking up the continuous 
mass and angular momentum loss into three discrete steps 
(see the left hand side of Fig. [l]; cf. also Langer 1998), ne- 
glecting secular changes of the stellar structure. Starting 
from a rigidly rotating envelope extending from R c to Rq 
(Al), the outer part of the convective envelope located 
between R\ and i?o is removed by stellar mass loss within 
the time interval Ai (Al — * A2). In the second step, for 
which we assume local angular momentum conservation, 
the envelope re-expands to roughly its original size (A2 — * 
A3). This results in a slow-down of the surface rotation 
and, as the layers deep inside the envelope expand less, 
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Fig. 1. Mass loss from a rigidly rotating stellar envelope from the surface (case A: left panels) and through its lower 
boundary (case B: right panels). For illustration, the continuous process is split up into three steps. First (panels 
1 — » 2), mass is removed from the envelope, second, the envelope restores its original radial extent through expansion 
(panels 2 — > 3), and third (panels 3 — > 4), the angular momentum is redistributed such that rigid rotation is restored. 
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lar momentum (A3 — > A4). As this implies an averaging 
of the angular velocity, it is clear from Fig 1 (A4) that 
now the whole convective envelope rotates slower than at 
the beginning of step 1. Obviously, the redistribution of 
angular momentum towards the stellar surface leads to an 
increase of the angular momentum loss rate. 

The efficiency of the angular momentum loss induced 
by mass loss from the surface of a rigidly rotating enve- 
lope, i.e. the amount of angular momentum lost per unit 
mass lost relative to average specific angular momentum 
of the envelope, is given by 

2 



J M„. 



X := 



Jsurf 



(r 2 



(2) 



where M cnv and J onv are total mass and total angular 
momentum of the envelope, M and J are stellar mass and 
angular momentum loss rate, j S urf and i S urf are the specific 
angular momentum and moment of inertia at the surface, 
respectively, and i?* is the stellar radius. The mean value 
of a quantity x over the envelope is defined by 

fM 



1 



Mr 



x(m) dm 

M-M onv 



where M is the mass of the star. The larger x the more 
efficient is the loss of angular momentum per unit mass 
lost. A value of x > 1 corresponds to a decrease of the 
mean specific angular momentum of the envelope, \ < 1 
to an increase. For the case of mass loss from the surface 
of a rigidly rotating stellar envelope x is always greater 
than 1. 

The density stratification in the envelope determines 
how much angular momentum is stored in the layers close 
to the surface relative to the total angular momentum of 
the envelope. An envelope structure which holds most of 
the mass close to its bottom favors high angular momen- 
tum loss rates, since this decreases (f 2 ) env - In red super- 
giants, on the other hand, (r 2 ) cnv is rather high, since 
those stars store a large fraction of their mass in layers far 
from the stellar center, thus reducing \. 

This is demonstrated in Fig. ||, where we plotted the 
moment of inertia per unit radius, dJ(r)/dr = 87rr 4 p/3 for 
spherical symmetry, as a function of the fractional radius 
for different types of stars. For a red supergiant model 
with extended convective envelope we find the major con- 
tribution to the total moment of inertia I of the star at 
radii around ~ 2/3 (x ~ 9/4), whereas for the zero-age 
main sequence (ZAMS) model matter at ~ 1/3 R+ domi- 
nates the moment of inertia of the star (x ~ 9) . As main 
sequence stars can be approximated by rigid rotators (cf. 
Zahn 1994), the whole star takes the role of the stellar 
envelope as far as our definition of x is concerned. There- 
fore, the considered main sequence star loses its angular 
momentum four times more efficient than the red super- 
giant model plotted in Fig. @. 




0.2 0.4 0.6 0.. 

radius / stellar radius 

Fig. 2. Moment of inertia per radius d/(r)/dr as a func- 
tion of radius (in units of the stellar radius) for four differ- 
ent stellar models. The solid line shows a ZAMS model, 
the dotted line a red supergiant model with an almost 
fully convective hydrogen-rich envelope, the dashed line a 
red supergiant model just before the blue loop, and the 
dash-dotted line a blue supergiant during the blue loop. 
All models are taken from the 12 M sequence described 
in detail below. 



supergiant stage, where I is dominated by the density in- 
version at the upper edge of the convective region, making 
the angular momentum loss from the envelope of this star 
quite inefficient (x ~ 1.5), and a blue supergiant model 
during central helium burning in which the moment of 
inertia is concentrated even more towards the center of 
the star (x » 30) than for the ZAMS model. Thus, if 
the envelopes of blue supergiants stay close to solid body 
rotation — which is indicated by our time-dependent cal- 
culations — they experience a more efficient spin-down 
than main sequence stars. 

However, since for blue supergiants the total moment 
of inertia / is smaller than for their progenitor main se- 
quen ce s tars, they may get closer to critical rotation (cf. 
Sect. |3.l| ) if they keep their angular momentum. This may 
be in particular the case for metal poor massive stars. 
Those have much smaller mass loss and therefore angu- 
lar momentum loss rates, and they can evolve from the 
main sequence directly into a long-lasting blue supergiant 
stage without an intermediate red supergiant phase (cf. 
e.g. Schaller et al. 1992); we found such stars to obtain 
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2.3. The spin-up of convective envelopes 

The spin-up of convective envelopes which decrease in 
mass due to an outflow through their lower boundary is 
understood in a similar way to the mass loss process dis- 
cussed above. Again we want to split up this process into 
tree discrete steps (see sequence B on the right hand side of 
Fig. 0). Starting from a rigidly rotating envelope extend- 
ing from Rq to R* (Bl), the inner part of the convective 
envelope, located between i?o and Ri, becomes radiative 
within some time At (B2). Then the envelope re-expands 
to its original size. As local angular momentum conser- 
vation is assumed in this step, it results in a spin-up of 
most of the convective envelope, which is weaker for the 
layers farther out. The outcome is a differentially rotat- 
ing envelope which spins fastest at its bottom (B3). Since 
the envelope has to go back to rigid rotation, angular mo- 
mentum is transported upward in the next step until this 
is achieved (B4). We end up with higher specific angular 
momentum in the whole convective envelope compared to 
the initial configuration. 

As can be seen in Fig. ^| below, during the evolution of 
a red supergiant towards the blue part of the HR diagram 
the radial extent of the convective envelope remains about 
fixed while mass shells drop out of the convective region. If 
one imagines the convective envelope to consist of moving 
mass elements or blobs, the spin-up process can thus also 
be understood as follows. A blob, starting somewhere in 
the convective region will, as it approaches the lower edge 
of the convective envelope, decrease its specific moment 
of inertia and therefore has to lose angular momentum in 
order to remain in solid body rotation with the whole con- 
vective region. Angular momentum has to be transferred 
to rising blobs such that also they remain in solid body 
rotation. Mass elements leaving the convective envelope 
thus only remove small amounts of angular momentum 
from the convective region. Therefore, the average specific 
angular momentum of the remaining convective envelope 
will increase and thus it spins up. 

Replacing j surf , z surf and R± in Eq. (§) by ji ow , i low and 
-Riow ! the specific angular momentum, the specific moment 
of inertia and the radius at the lower edge of the convective 
envelope, respectively, one can define an efficiency of an- 
gular momentum loss through the lower boundary, \- We 
find x "C 1, especially for those cases where x is small, as 
e.g. for the red supergiant envelopes under consideration 
here. 

The total change of specific angular momentum of the 
envelope (j) onv by mass loss through the upper and lower 
boundary can the be written as 

din Q') cnv = Af(l-x)+A/iow(l-x) (3) 
dt Af cnv 

where M\ ow is the rate at which mass leaves the envelope 



supergiant stage preceding the evolution off the Hayashi 
line, so that (j) almost increases as M env decreases, 

dln 0')cnv/ d lllM onv > - I , 

reflecting the fact that angular momentum does not get 
lost efficiently from the convective envelope through its 
upper nor its lower boundary. The approximation x« 1 
is hardly affected by the variation of the lower bondary 
radius seen in Fig. ||. Thus, the convective envelope loses 
most of its mass but keeps a major part of the angular 
momentum. 

When the convective envelope gets depleted in mass 
and the stellar radius decreases considerably, the global 
contraction of the stellar envelope results in an additional 
contribution to its spin-up. Still, mass elements drop out 
of the convective region (cf. Fig. ||), but now the con- 
traction leads to an increase of the rotation velocity, and 
the star can reach rotation velocities of the order of the 
break-up velocity (see below). However, the contraction 
does not contribute to the increase of the specific angular 
momentum at the surface. 

Finally, we want to note that if the whole star would 
remain rigidly rotating, e.g., due to the action of magnetic 
fields inside the star or by more efficient shear instabilities, 
its spin-up would occur very similar to the case described 
here. In that case, also mass shells from the core would 
transfer part of their angular momentum to layers above, 
which would make the spin-up somewhat more efficient. 
However, since in a red supergiant the mass elements lose 
the major part of their angular momentum to the con- 
vective envelope before they leave it, and the amount of 
angular momentum contained in the core is small anyway 
(cf. Sect. the additional spin- up will be small (cf. 

Fig. before blue loop). 

3. Numerical simulations 

3.1. Input physics 

The stellar evolution calculations presented here are ob- 
tained with an implicit hydrodynamic stellar evolution 
code (cf. Langer et a!. 1988). Convection according to the 
Ledoux criterion and semiconvection are treated accord- 
ing to Langer et al. (1983), using a mixing length param- 
eter «mlt = 1.5 and semiconvective mixing parameter of 
d sem = 0.04 (Langer 1991a). Opacities are taken from 
Alexander (1994) for the low temperature regime, and 
form Iglesias et al. (1996) for higher temperatures. The 
effects of rotation on the stellar structure as well as the 
rotationally induced mixing processes are included as in 
Pinsonneault et al. (1989), with uncertain mixing efficien- 
cies adjusted so as to obtain a slight chemical enrichment 
at the surface of massive main sequence stars (cf. Fliegner 
et al. 1996). 
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tationally induced enhancement of mass loss as the star 
approaches the il-Limit (cf. Langer 1997), i.e. 



M(u>) 



1 



M(w = 0) V 1 -^, 
(cf. Friend & Abbott 1986) where 



n := 



Gi 



(i-r) 



r := 



incGM 



Here, k is the Rosseland opacity. T is not only evaluated at 
the stellar surface, but its maximum value in the radiative 
part of the optical depth range r £ [2/3, 100] (cf. Lamers 
1993; Langer 1997) is used. 

The quantitative result for the ^-dependence of the 
mass loss rate obtained by Friend & Abbott (1986) was 
questioned by Owocki et al. (1996), who performed hydro- 
dynamic simulations of winds of rotating hot stars includ- 
ing the effect of non-radial radiation forces and gravity- 
darkening in the way described by von Zeipel (1924). How- 
ever, the only crucial ingredient for the model calculations, 
which is confirmed by Owocki et al. (1997), is the fact that 
the mass loss rate increases strongly as the star approaches 
the n-limit, so that the star cannot exceed critical rota- 
tion, but rather loses more mass and angular momentum 
(cf. also Langer 1998). 

To quantify the angular momentum loss due to stellar 
winds, J, we assume that, at any given time, the mean 
specific angular momentum of the wind material leaving 
the star equals the specific angular momentum averaged 
over the rigidly rotating, spherical stellar surface, which 
we designate as j sur f . 

The transport of angular momentum inside the star is 
modeled as a diffusive process according to 



du> 
~dt 
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d 



Airr 2 pi dr 



a 2 • n ® U 

4irr piD—- 
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where D is the diffusion coefficient for angular momen- 
tum transport due to convection and rotationally induced 
instabilities (cf. Endal & Sofia 1979; Pinsonncault et al. 
1991), and the last term on the right hand side accounts 
for advection. The diffusion equation is solved for the 
whole stellar interior. In stable layers, the diffusion co- 
efficient is zero. We specify boundary conditions at the 
stellar surface and the stellar center which guarantee an- 
gular momentum conservation to numerical precision. 

Our prescription of angular momentum transport 
yields rigid rotation when and wherever the time scale 
of angular momentum transport is short compared to the 
local evolution time scale, no matter whether rotationally 
induced mixing or convective mixing processes are respon- 
sible for the transport. Chaboyer & Zahn (1992) found 
that meridional circulations may lead to advection terms 
in the angular momentum transport equation, which may 



main sequence phase (cf. Talon et al. 1997). However, 
the omission of these terms has no consequences for the 
spin-up process describe here, since it is dominated by the 
much faster angular momentum transport due to convec- 
tion. Since the time scale of convection is generally much 
smaller than the evolution time, convective regions are 
mostly rigidly rotating in our models. 

The effect of instabilities other than convection on the 
transport of matter and angular momentum are of no rel- 
evance for the conclusions obtained in the present paper. 
E.g., mean molecular weight gradients have no effect on 
the spin-up process described below since it occurs in a 
retreating convection zone (see below). 

Calculations performed with a version of the stellar 
evolution code KEPLER (Weaver et al. 1978), which was 
modified to include angular momentum, confirm the spin- 
up effect obtained with our code which is described in 



Sect. 3.2 



3.2. Results 

We have computed stellar model sequences for different 
initial masses and compositions (cf. Heger et al. 1997), 
but here we focus on the calculation of a 12 M© star of 
solar metallicity. This simulation is started on the pre- 
main sequence with a fully convective, rigidly rotating, 
chemically homogeneous model consisting of 28 % helium 
and 70 % hydrogen by mass and a distribution of metals 
with relative ratios according to Grevesse & Noels (1993). 
Its initial angular momentum of 110 1 50 ergs leads to an 
equatorial rotation velocity of ~ 200 kms -1 on the main 
sequence, which is typical for these stars (cf. Fukuda 1982; 
Halbedel 1996). 

During the main sequence evolution the star loses 25 % 
of its initial angular momentum and 0.23 M of its enve- 
lope due to stellar winds. After the end of central hydrogen 
burning, the star settles on the red supergiant branch af- 
ter several 10 yr and develops a convective envelope of 
8.3 M a . During this phase it experiences noticeable mass 
loss (M ss 8 1O _7 M0 yr -1 ), but the angular momentum 
loss per unit mass lost J/M w 1.7 10 19 cm 2 s _1 is lower 



than on the main sequence (cf. Sect. 2.2). 

Shortly after core helium ignition, the bottom of the 
convective envelope starts to slowly move up in mass. 
About 4 10 5 yr later (at a central helium mass fraction 
of 65%), the convective envelope mass decreases more 
rapidly. After another ~ 5 10 5 yr it reaches a value of 
~ 3M Q . Up to this time, the star has lost 0.75 M Q as 
a red supergiant, and 43 % of the angular momentum left 
at the end of the main sequence. 

At its largest extent in mass, the convective envelope 
contains a total angular momentum of 77 1 50 ergs, of 
which 17 10 50 ergs are contained in the lower 4.5 M Q . Af- 
ter those layers dropped out of the convective envelope, 
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Fig. 3. Track of the 12 M Q star in the HR diagram during 
its blue loop. The two arrows indicate the direction of the 
evolution. The thick drawn part of the track corresponds 
to the time span shown in Figs. || and |?]. The star spends 

1 000 yr between two neighboring tick marks. The thick 
dot indicates where the star obtains its maximum mass 
and angular momentum loss rates; see the dotted line in 
Fig. |. 

lost due to mass loss. Thus, on average the specific angu- 
lar momentum of the lower part of the convective envelope 
is decreased by up to a factor of ~ 5, while that of the 
remaining convective envelope increased on average by a 
factor of 1.5, despite the angular momentum loss from the 
surface (cf. Fig. g Eq. @). For comparison, the helium 
core contains never more than ~ 1.5 10 50 ergs. 

At this point of evolution, the transition to the blue 
sets in in our model. Within the next ~ 25 000 yr another 

2 Mq drop out of the convective envelope which then com- 
prises only about 1 M Q but has still the full radial extent of 
a red supergiant. The ensuing evolution from the Hayashi 
line to the blue supergiant stage takes about 10 000 yr (cf. 
Fig. [|). During this time, the angular momentum trans- 
port to the outermost layers of the star continues, tapping 
the angular momentum of those layers which drop out of 
the convective envelope. 

The contraction of the star by a factor of / ss 10 (cf. 
Fig. ||) reduces the specific moment of inertia at the sur- 
face by f 2 » 100 and u would be increased by the same 
factor if j were conserved locally. This would imply an in- 
crease of the equatorial rotational velocity by a factor of 
- 10 (cf. Fig. fR "decoupled"). The true increase in the 
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Fig. 4. Evolution of the internal specific angular momen- 
tum profile in the hydrogen-rich envelope of the 12 M© se- 
quence during the first part of core helium burning which 
is spent as a red supergiant. Dashed lines show specific an- 
gular momentum versus mass coordinate in the convective 
part of the envelope for different times, from 420 000 yr 
to 2 000yr before the red — > blue transition, with those 
reaching down to lower mass coordinates corresponding 
to earlier times and extending to higher mass coordinates, 
because the stellar mass decreases with time. The solid 
line shows the angular momentum profile for those mass 
zones which have dropped out of the convective envelope, 
which remains frozen in later on. The dotted line marks 
the angular momentum profile in the inner stellar region 
which is never part of the convective envelope. 

be even larger if the star would not arrive at the f2-limit 
(see Fig. and lose mass and angular momentum at an 
enhanced rate. 

For local angular momentum conservation Q scales as 
n cx l/y/R(l - T). A value of T » 2 10~ 3 is found on the 
red supergiant branch in our calculation, and T « 0.4 for 
the blue supergiant stage. This by itself would increase 

by a factor of ~ 3.9 during the red — > blue transition. 
Actually, £1 changes from w 0.01 at the red supergiant 
branch at the beginning of central helium burning to Q sa 

1 — to critical rotation — at the red — > blue transition, 
i.e. by a factor ~ 100, despite significant mass and angular 
momentum loss. 

At an effective temperature of ~ 6 000 K the Edding- 
ton-factor T at the surface (as defined above) rises from a 
few times 10~ 3 to ~ 0.4, mainly due to an increase in the 
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500 r 




Fig. 5. Evolution of the radii of different mass shells as 
a function of time for a period including the transition 
from the red to the blue supergiant stage of the 12 M© 
model. The zero-point on the x-axis corresponds to the 
thick dot in Fig. || and corresponds roughly to the time 
of the red — ► blue transition. Except for the uppermost 
solid line, which corresponds to the surface of the star, 
the lines trace Lagrangian mass coordinates. The mass 
difference between the lines is 0.5 M Q . Shading indicates 
convective regions. 

centrated in the upper 0.01 M Q . Since becomes close 
to 1, the mass loss rate rises to values as high as several 
10~ 5 M Q yr -1 (cf. Fig. ||) in order to keep the star below 
critical rotation, with the result that these layer are lost 
within a few 1 000 yr. This is by far the most dramatic 
loss of angular momentum, i.e. the highest value of J, the 
star ever experiences. The specific angular momentum loss 
j/M reaches a peak value of 7.7 10 19 cm 2 s" 1 (cf. Fig. §). 

After arriving at the blue supergiant stage, the star 
still experiences a high angular momentum loss rate for 
some time, since the major part of the star's angular mo- 
mentum is still concentrated in the vicinity of the surface. 
Interestingly, the star now even spins faster than when 
rigid rotation of the whole star were assumed (cf. Fig. fjj), 
because within the blue supergiant 's radiative envelope 
the angular momentum is not transported downward effi- 
ciently. However, due to mass loss, this deviation does not 
persist long. 

Due to the long duration of the blue supergiant phase, 
several 10 5 yr in comparison to the 10 4 yr the red — > blue 
transition takes, the total angular momentum J is reduced 
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Fig. 6. Evolution of characteristic stellar properties as 
function of time, during the first part of the blue loop 
of our 12 M model (cf. Fig. |^). The time zero-point is 
defined as in Fig. || and is marked by the dotted line. 
Displayed are: the angular momentum loss rate J (A), 
the total angular momentum J (B), the mass loss rate M 
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Fig. 7. Equatorial rotation velocity as function of time 
(solid line) compared to the Keplerian (dashed line) and 
the critical (dotted line) rotation rate; the latter two are 
different by the factor y/l — T. During the red supergiant 
phase it is r <C 1 and the two lines coincide, while during 
the blue supergiant phase T rises to 0.4. The dash-dotted 
line shows the evolution of the surface rotation rate if there 
were no angular momentum transport in the convective 
envelope, and the dash-triple-dotted line shows how the 
surface rotation rate would evolve if the whole star would 
maintain rigid rotation. 

In total, after the blue loop the star has ~ 5 times less 
angular momentum than before, and thus red supergiants 
which underwent a blue loop rotate significantly slower 
than those which did not. 



it could initially be far away from critical rotation. Due to 
the spin-up mechanism described here, such a star might 
then approach critical rotation even before the transition 
into a blue supergiant. 

The mechanism to spin- up stellar envelopes (Sect. ^|), 
which was applied to a post-red supergiant 12 M Q star 
in Sect. g], may operate in all evolutionary phases dur- 
ing which stars evolve from the Hayashi-line towards the 
hotter part of the HR diagram. These comprise the transi- 
tion of pre-main sequence stars form their fully convective 
stage to the main sequence, the transition of low and in- 
termediate mass stars from the Asymptotic Giant Branch 
(AGB) to central stars of planetary nebulae, the blue loops 
of stars in the initial mass range ~ 3M Q ... ~ 25 M©, and 
the transition of massive red supergiants into Wolf-Rayet 
(WR) stars. For all four phases, observational evidence for 
axisymmetric circumstellar matter exists. 

A high value of f2 during the evolution off the Hayashi 
line can be expected to have notable influence not only on 
the mass and angular momentum loss rate as in the case 
studied in Sect. ||, but also on the geometry of the wind 
flow. At present, theoretical predictions of the latitudinal 
dependence of stellar wind properties for rapidly rotating 
stars are ambiguous. Bjorkman & Cassinelli (1993) have 
proposed that angular momentum conservation of parti- 
cles in a pressure free wind which is driven by purely radial 
forces leads, for sufficiently large values of O, to high ratios 
of equatorial to polar wind density, i.e. to disks or disk- 
like configurations. Owocki et al. (1996) found that this 
might no longer be the case when the non-radial forces 
occurring in hot star winds are accounted for. For cooler 
stars, the prospects of disk formation may thus be better 
(cf. also Ignace et al., 1996). To discus the geometry of 
the winds of our models is beyond the scope of this paper; 
however, we want to point out that the maximum values 
of Q and j/M occur at about T cS ~ 6 000K (cf. Sect. |1 
and Fig. ||) , so that an equator-to-pole wind density ratio 
larger than one may still be justified. 



4. Discussion 

The mechanism to spin-up stellar envelopes presented 
above enhances the specific angular momentum in the 
surface layers of the star by some factor (cf. Sect. |3.2| ); 
a question to be considered separately is the origin of the 
angular momentum, and also in which cases critical ro- 
tation is reached. For single stars, the available angular 
momentum seems to be limited to that of the ZAMS star, 
reduced by the angular momentum lost during its evolu- 
tion. A way to supply a considerable amount of angular 
momentum to the star would be the capture of a com- 
panion star or planet (Soker 1996; Podsiadlowski 1998). 
This might occur when the star first becomes a red su- 
pergiant. The resulting configuration would have a much 



4--1- Blue loops during or after core helium burning 

An interesting example where the spin-up mechanism de- 
scribed in this paper should almost certainly have played 
a role is the progenitor of SN 1987A. In fact, the SN 1987A 
progenitor is the only star of which we are reasonably sure 
that it performed a blue loop, in this case after core helium 
exhaustion (cf. Arnett et al. 1989) 

The structures observed around SN 1987A seem to be 
rotationally symmetric with a common symmetry axis, 
which may suggest that rotation has played a major role 
in their formation. The inner of the three rings is cur- 
rently explained by the interaction of the blue supergiant 
wind with the wind of the red supergiant precursor (cf., 
Chevalier 1996). However, this interaction would result in 
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gion, and maybe also the two outer rings, significant rota- 
tion appears therefore to be required (cf. Martin & Arnett 
1995), which may be provided by our spin-up mechanism. 
We may note that for this mechanism to work it is insignif- 
icant what actually triggered the blue loop; i.e. it would 
work as well for single star scenarios with a final blue loop 
(cf. Langer 1991b; Meyer 1997; Woosley et al. 1998) or for 
binary merger scenarios which predict a final red — > blue 
transition of the merger star (cf. Podsiadlowski 1998). 

We do not expect to find ring nebulae frequently 
around blue supergiants (cf. Brandner et al. 1997), since 
the time scale on which our model (cf. Sect. ||) shows very 
high surface rotation rates is rather small compared to the 
typical life time of a blue supergiant (cf. Langer 1991a), 
and because they are quickly dissolved by the blue su- 
pergiant wind. SN 1987A is an exception, since here the 
transition to the blue happened only short time ago (i.e. 
the supernova exploded only shortly after the transition). 
However, a certain type of B supergiants, the B[e] stars, 
show emission line features which might be due to a cir- 
cumstellar disk (cf. Zickgraf et al. 1996). The location of 
the less luminous subgroup of B[e] stars in the HR diagram 
(Gummersbach et al. 1995) is in fact consistent with a blue 
loop scenario for their evolution. I.e., their disks might be 
produced by the spin-up mechanism described here (cf. 
also Langer & Heger 1998). 



4-2. Red supergiant — > Wolf-Rayet star transition 

The transition of massive mass losing red supergiants into 
Wolf-Rayet stars is the massive star analogue of the AGB 
— > post-AGB star transition. It occurs when the mass 
of the hydrogen-rich envelope is reduced below a critical 
value (cf., e.g., Schaller et al. 1992) and should not be con- 
fused with the blue loops discussed in Sect. LI. As in the 
post-AGB case (Sect. L4), the spin- up mechanism ca n be 
expected to work. Unlike in the case discussed in Sect. 3.2, 
the major part of the envelope is lost due to stellar winds 
before the red supergiant — > Wolf-Rayet star transition, 
with the consequence of considerable angular momentum 
loss. Therefore, it may be more difficult for the star to 
reach critical rotation during the contraction. 



4-3. Pre-main sequence evolution 



As mentioned in Sect. 3.2, we started our 12 M Q sequence 



However, there are signs of asphericity in the ring 
nebula NGC 6888 around the Galactic Wolf-Rayet star 
WR 136 which has been interpreted as swept-up red super- 
giant wind shell by Garci'a-Segura & Mac Low (1995) and 
Garcfa-Segura et al. (1996). Also, Oudmaijer et al. (1994) 
report on bi-polar outflows from IRC+10420, a massive 
star just undergoing the red-supergiant — * Wolf-Rayet star 
transition. IRC+10420 is currently an F type star, i.e. it 
has an effective temperature at which we expect the max- 
imum effect of the spin-up mechanism discussed here (cf. 
Sect. ETl). Therefore, in the absence of a binary compan- 



from a fully convective pre-main sequence configuration. 
We expected the spin-up mechanism found for the post- 
red supergiant stage to be also present in the contrac- 
tion phase towards the main sequence. An analysis of this 
evolutionary phase showed in fact its presence, although 
the efficiency of the spin-up was found to be somewhat 
smaller. During that part of the pre-main sequence con- 
traction phase where the convective region retreats from 
the center of the star to its surface, the star would have 
increased its equatorial rotational velocity by a factor of 
~ 3 if angular momentum were conserved locally, but it 
was spun- up by a factor of ~ 10. I.e., the spin- up mecha- 
nism described in Sect. || resulted in an additional increase 
of the rotational velocity of a factor of ~ 3. This was 
not enough to bring the star to critical rotation in this 
phase. However, if the initial rotation rate would have 
been larger, a phase of critical rotation during the pre- 
main sequence stage would well be possible. 

While a fully convective pre-main sequence stage may 
not be realistic for massive stars and is just assumed 
in our case for mathematical convenience, pre-main se- 
quence stars of low and intermediate mass are supposed 
to evolve through a fully convective stage (Palla & Stahlcr 
1991). During the transition from this stage to the main 
sequence, the spin-up mechanism described in this paper 
might operate. Pre-main sequence stars in the correspond- 
ing phase, i.e. past the fully convective stage but prior 
to core hydrogen ignition, are often found to have disks. 
They correspond to the T Tauri stars at low mass (e.g., 
Koerner & Sargent 1995) and to central stars of Herbig 
Haro objects at intermediate mass (e.g., Marti et al. 1993). 
However, the disks are usually interpreted as remnants of 
the accretion process which built up the star. Since pre- 
main sequence stars are often found to be rapid rotators 
(cf. Walker 1990), we may speculate here about a possible 
contribution to the disk due to decretion from the star 
reaching critical rotation due to spin-up (cf. Krishnamur- 
tihi et al. 1997). 

4--4- Post-AGB evolution 

Low and intermediate mass stars leave the AGB when the 
mass of their hydrogen-rich envelope decreases below a 
certain value. When this happens, the envelope deflates, 
the stellar radius decreases, and the energy transport in 
the envelope changes from convective to radiative. The 
spin-up mechanism described in Sect. can be expected 
to operate in this situation, as in the red supergiant — > 
Wolf-Rayet star transition (Sect. p~2| ). 

Whether or not post-AGB stars can reach critical ro- 
tation due to this spin-up is not clear. Certainly, the heavy 
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The ratio of the rotation rate to the critical rotation rate f2 
may be further affected by random kicks due to asymmet- 
ric mass loss, which may keep the envelope at some level 
of rotation (cf. Spruit 1998), or by transport of angular 
momentum from the core to the envelope, which may be 
efficient during the thermal pulses, and by the evolution of 
the critical rotation rate, which depends on the Edding- 
ton factor V and thereby on the opacity coefficient (cf. 
Garcfa-Segura et al. 1998). 

Clearly, the spin-up mechanism described in this paper 
may help to bring post-AGB stars — or more specific: 
stars which just left the AGB, i.e. central stars of proto- 
planetary nebulae — closer to critical rotation. It may play 
a role in explaining axisymmctric flows which are often 
observed in central stars of proto-planetary nebulae (cf. 
Kwok 1993) and the shapes of bi-polar planetary nebulae 
(cf. Schwarz et al. 1993; Stanghellini et al. 1993; Garcfa- 
Segura et al. 1998). 

5. Summary and conclusions 

In this paper, we discussed the effect of mass outflow 
through the inner or outer boundary of a rigidly rotat- 
ing envelope on its rotation frequency. It causes a change 
of the specific angular momentum in the envelope and al- 
ters its rotation rate besides what results from contraction 
or expansion (cf. Fig. 0). For constant upper and lower 
boundaries of the envelope, which we found a good ap- 
proximation for convective envelopes (cf. Fig. |5|), a spin- 
down occurs for mass outflow through the upper boundary 
— which corresponds, e.g., to the case of stellar wind mass 
loss from a convective envelope (cf. also Langer 1998) — , 
while a spin-up results from mass outflow through the 
lower boundary (cf. Fig. |l|). 

The latter situation is found in evolutionary models of 
a rotating 12 Mq star at the transition from the Hayashi- 
line to the blue supergiant stage. The star increased its ro- 
tational velocity one order of magnitude above the velocity 
which would result in the case of local angular momentum 
conservation. It would have increased its rotational veloc- 
ity even further if it would not have arrived at critical 
rotation (cf. Sect. IL2, Fig. 0), with the consequence of 
strong mass and angular momentum loss. At this point, 
the specific angular momentum loss J/M reached about 
810 19 cm 2 s- 1 (cf. Fig. |). 

The geometry of circumstellar matter around stars 
which undergo a red — > blue transition may be strongly af- 
fected by the spin-up. We propose that this was the case 
for the progenitor of SN 1987A, the only star of which 
we know that it performed a red — > blue transition in the 
recent past. The blue supergiant in its neighborhood stud- 
ied by Brandner et al. (1997), around which they found a 
ring nebula as well, is another candidate. Also, B[e] stars 
may be related with the post-red supergiant spin-up (cf. 



central stars of proto-planetary nebulae (Sect. 4.4), from 
stars in the transition phase from the red supergiant stage 
to the Wolf-Rayet s tage (Sect. 4.2), and from pre-main se- 
quence stars (Sect. 4.3). 
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